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We compute the strangeness and light-quark contributions As, Au and Ad to the proton spin 
in 7if = 2 lattice QCD at a pion mass of about 285 MeV and at a lattice spacing a ~ 0.073 fm, 
using the non-perturbatively improved Sheikholeslami-Wohlert Wilson action. We carry out the 
renormalization of these matrix elements which involves mixing between contributions from different 
quark flavours. Our main result is the small negative value As MS {V7A GeV) = -0.020(10)(4) of 
the strangeness contribution to the nucleon spin. The second error is an estimate of the uncertainty, 
due to the missing extrapolation to the physical point. 
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Introduction. — The proton spin can be split into a 
quark spin contribution AS, a quark angular momentum 
contribution L q and a gluonic contribution AG (includ- 
ing spin and angular momentum) [lj: 



\ = \^ + L q + AG. 



(1) 



In the naive non-relativistic SU(6) quark model, AS = 1, 
with vanishing L q and AG. In this case there will be no 
strangeness contribution As to AE = Au+Ad+As+- ■ ■ , 
where, in our notation, Aq — AE g contains both, the 
spin of the quarks q and of the antiquarks q. 

Experimentally, As is obtained by integrating the 
strangeness contribution As(x) to the spin structure 
function g\ over the momentum fraction x. The inte- 
gral over the range in which data exist agrees with zero; 
see, e.g., new COMPASS data @, Q for x > 0.004 or 
HERMES data 3 for x > 0.02, while global analyses 
give values [H|— [3] As ~ —0.12, suggesting a large negative 
As(x) at very small x. Pioneering lattice simulations of 
disconnected matrix elements also indicated values 0, 0] 
As w —0.12. However, the errors given in these studies 
are quite optimistic while the global fits rely on an ex- 
trapolation of the integrated experimental AE to small 
x and constrain the axial octet charge as to a value, ob- 
tained from hyperon /3-decays, assuming SU(3) f flavour 
symmetry. Some time ago, employing heavy baryon chi- 
ral perturbation theory, Savage and Walden [Io[ pointed 
out that SU(3)f symmetry in weak baryonic decays may 
be violated by as much as 25 % and hence As (a;) could 
remain close to zero also for x < 0.001; see also [il| . 
SU(3) f symmetry is however supported by lattice simu- 



lations of hyperon axial couplings 12 - lB[ , albeit within 
non-negligible errors. 

In this Letter, we directly compute the matrix elements 
that contribute to the Ag, including quark line discon- 
nected diagrams. Preliminary results were presented at 
conferences 

Simulation details and methods. — We simulate nt = 
2 non-perturbatively improved Sheikholeslami-Wohlert 
Fermions, using the Wilson gauge action, at j3 — 5.29 
and k = k u( i = 0.13632. Setting the scale from the chi- 
rally extrapolated nucleon mass [l9j], we obtain the lat- 
tice spacing a -1 = 2. 71 (2) (7) GeV, where the errors are 
statistical and from the extrapolation, respectively. 

We realize two additional valence n values, K m = 
0.13609 and n s — 0.13550. The corresponding 
pion masses are mps.ud = 285(3)(7) MeV, mps, m = 
449(3)(ll)MeV and m PStS = 720(5)(18) MeV. k s was 
fixed so that the mps, s value is close to the mass of 
a hypothetical strange-antistrange pseudoscalar meson: 



(m 2 K ± 



ml ± ) 



1/2 



686.9 MeV. We investi- 



gate volumes of 32 J 64 and 40 J 64 lattice points, i.e., 
Lm-ps.ud — 3.36 and 4.20, respectively, where the largest 
spatial lattice extent is L s» 2.91 fm. 



The quark polarizations are extracted from the large- 
time behaviour of ratios of three-point over two-point 
functions. We create a polarized proton at a time to = 0, 
probe it with an axial current at a time t and destroy 
the zero momentum proton at tf > t > 0. Quark line 
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connected and disconnected terms contribute: 



R con (tut) 



R dis (t u t) 



pol Spt 



(*)> 



(2) 



ClO*f) Ex TrbnsM-^x, t; x, t)]> 



/pa/3 

\ unpol 2pi 



Ct(tt)) 



Here M is the lattice Dirac operator, r unpo i = |(1 + 74) 
is a parity projector and r po i = 27j75r unpo i projects out 
the difference between the two polarizations (in direction 
3). We average over j = 1, 2, 3 to increase statistics. For 
the up and down quark matrix elements we compute the 
sum of connected and disconnected terms while only R dls 
contributes to As. 

For disconnected contributions we fix the time distance 
between the source and the current insertion t = 4a ~ 
0.29 fm and vary tf. Both t and the distance between 
current and sink tf — t should be taken large, to suppress 
excited state contributions. Using the sink and source 
smearing described in [20j . we find the asymptotic limit 
to be effectively reached for tf ~ 6a-7a; see Fig. [T] for an 
example. The saturation into a plateau at tf < 2t and the 
convergence of the point sink data towards the same value 
demonstrate that t — 4a was reasonably chosen. To be on 
the safe side, we only fit the tf > 8a f« 0.58 fm smeared- 



smeared ratios. Building upon previous experience 21|, 
the connected part, for which the statistical accuracy is 
less of an issue, is obtained at the larger, fixed value 
tf = 15a, varying t. 

The disconnected contribution is computed with the 
stochastic estimator methods described in em- 
ploying time partitioning, a second order hopping pa- 
rameter expansion and the truncated solver method. We 
compute the Green functions for four equidistant source 
times on each gauge configuration. We also construct 
backwardly propagating nucleons, replacing the positive 
parity projector ^(1 + 74) by i(l — 74), seeding the 
noise vectors on eight (four times two) timeslices. In 
addition to the 48 (four times spin times colour) solves 
for smeared conventional sources, that are necessary to 
construct the two-point functions, we run the Conjugate 
Gradient (CG) algorithm on Ni = 730 complex 1*2 noise 
sources for n t = 40 iterations. The bias from this trunca- 
tion is corrected for 22j by N2 = 50 BiCGstab solves that 
are run to convergence. We analyse a total of 2024 ther- 
malized trajectories on each of the two volumes where we 
bin the data to eliminate autocorrelations. 

Renormalization. — Non-singlet axial currents renor- 
malize with a renormalization factor Z r A {a) that only 
depends on the lattice spacing. This was determined 
non-perturbatively for the action and lattice spacing in 
use = 0.76485(64)(73). 

However, due to the axial anomaly, the renormaliza- 
tion constant of singlet currents, Z A (/j,,a), acquires an 
anomalous dimension. To first non-trivial order this 
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FIG. 1. The disconnected ratio R 
volume at « va i = «: cur = k s for smeared-smeared (SS) 
smeared-point (SP) source-sink combinations. 



and 



reads |2J, 125J j s A {a s ) = -6C F n f [a s /(47r)] z 



Z S A deviates 

from Z™ starting at O(oj.) in perturbation theory. Both 
factors have been calculated to this order, with the result 
for the conversion into the MS scheme at a scale /i 



2(5] 



z( f x,a) = Z A ( f i 1 a)-Z A s (a) 

= C F n f [15.8380(8) - 61n(aV)] {^j 2 , (3) 

where we have set the Sheikholeslami-Wohlert parame- 
ter csw — 1 to be consistent to this order in perturba- 
tion theory. To this first non-trivial order no scale enters 
the coupling parameter a s . Since perturbation theory 
in terms of the bare lattice parameter ao — 6/ (47T/3) 
is known to converge poorly, we substitute a s by a 
coupling defined from the measured average plaquette 
a s = -31n(J7 n )/(47r) = 0.14278(5), where we have used 
the chirally extrapolated value [23 (U n ) = 0.54988(11). 

No dimension-four operator can be constructed that 
mixes with the relevant forward matrix element of 97^75 q 
and that cannot be removed, using the equations of mo- 
tion [2^. This also holds for the singlet case 29], such 
that we only need to replace 



i->- 



Z n A s {l + b A am), Z\^t Z s A (l + b A am), (4) 



to achieve full 0(a) improvement. The factor b A is known 
to 0(a s ) [3: b A = b\ + 0{a 2 s ) w 1 



We obtain the values 



18.02539 C F ^. 

r 47T 



1 + b A am 



1.0324(3)(47) 
1.0041(3)(5) 



(m s , k 
(m ud , n 



0.13550) 
0.13632) 



(5) 



where the first error is due to the uncertainty in the quark 
mass and the second error corresponds to 50% of the 
one-loop correction. Considering the small size of this 
correction it is unlikely that the (two-loop) difference be- 
tween singlet and non-singlet ^-factors will result in any 
noticeable effect, and in particular not at the light-quark 
mass m u d, where it will be needed [see Eq. (|lll) below]. 
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ror to allow for higher order corrections. Due to the 
small anomalous dimension that only sets in at Q(a s ), 
the difference between singlet and non-singlet renormal- 
ization constants remains small, also at other scales. 
For instance, we obtain z(VlOGeV) = 0.0049(25) and 
z(2GeV) = 0.0082(41). 

In the nf = 1 + 1 + 1 theory the matrix elements renor- 
malize as follows: 



(m ps ) (val)(GeV) 

FIG. 2. Volume and valence quark mass dependence of the 
unrenormalized As lat . 



For m = 2 we get 

z{VTaGc\) = 0.0055(1)(27) , (6) 

at the renormalization scale p? = 7.4 GeV 2 = 
1.01(5) a -2 . We again include a 50 % systematic er- 

I 



g A = AT 3 = (Au - Ad) 



MS 



= Zf(a)(Au- Ad) lat (a) 



AT 8 = (Au + Ad - 2As) 



MS 



Z^(a)(Au + Ad -2 As 



,lat/ 



a = AE MS (u) = (Au + Ad + As) MS (/i) 



MS/ 



(7) 



(8) 



= Z s A (ii,a)(Au + Ad + A S y a \a). (9) 

We remark that for non-equal quark masses the non- 
singlet combinations Eqs. Q and © also receive con- 
tributions from disconnected quark line diagrams. 

We employ m = 2 sea quarks so that our singlet cur- 
rent is Au + Ad rather than the AS of Eq. ©. This 
modifies the renormalization pattern: 



MS 



Au(/x) \ 
Ad(fj) 
As(/x) / 



z(jU,q) 



z(/i,a) 
2 

2 



z(/4,a) 
2 



Zl«(a) 



2 







ZJ(a) 



/ Au(a) 
Ad(a) 
\ As (a) 



lat 



(10) 



As receives light-quark contributions but the Au MS 
and Ad MS remain unaffected by the (quenched) strange 
quark. Obviously, unitarity is violated, due to this 
quenching. The combination AT% still transforms with 
Z% s [Eq. ©] while Eq. © is violated, as it should be; 
instead, the Uf = 2 singlet operator Au + Ad renormal- 
izes with Z S A . We remark that the above renormalization 
pattern is similar to that of the scalar matrix element in 
the nf = 2 theory [HEmi. Note that in spite of the 
quenched strange quark the mismatch between directly 
converting the result into the MS scheme at a scale fi, us- 
ing z{nf = 2)/2, and first converting into the MS scheme 
at another scale // and subsequently running within the 
MS scheme with ln(/i///)7^(rif = 3)/3 to the scale fi is 
tiny. 

Results and systematics. — In Fig. [5] we display the vol- 
ume and (light) valence quark mass dependence of our 
unrenormalized As lat . There are no statistically signifi- 
cant finite size or mass effects. 

Using Eqs. (flU|) and ((H) we can renormalize 



Ag Mb (/i) =Z^(l + b A am q )Aq 



lat 



(Au + Ad) 



lat 



(11) 



for q G {u, d, s}. As discussed above, we omit the 
0(a) improvement factor (b s A Z A — b A Z I A )am u d of the 
(Au + A<i) lat term. This is of 0(a 2 s am u d) and numer- 
ically negligible. We display the bare lattice numbers 
for the connected and disconnected contributions to the 
proton spin and the renormalized 0(a) improved values 
in Table [fl for the two volumes. The Au^ and Ad m 
values are reduced by about 0.035, due to the sea quark 
contributions while As MS increases by 0.002 (< 10 %), 
due to the mixing with light-quark flavours. 

The uncertainties associated to the renormalization 
are much smaller than the statistical errors. Below 
we will only quote large volume results, with statistical 
and renormalization errors added in quadrature. Error 
sources that have so far not been accounted for are the 
missing continuum limit extrapolation, the quenching of 
the strange quark and simulating at a light sea quark 
mass value that is four times bigger than the physical 
one. There are no indications of radical quark mass ef- 
fects: the flavour mixing effects within the renormaliza- 
tion are small in spite of the comparatively large Au and 
Ad values. The dependence on the valence quark mass 
is small too; see Fig. [5J 
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TABLE I. The connected and disconnected contributions to 
Aq lat as well as the renormalized spin content at a scale /i = 
\/7.4 GeV. (The ATi are scale- independent.) The first error is 
statistical, the second is from the renormalization. In addition 
an overall 20 % systematic error needs to be added. 



q 




V,L 


^ycon 




Ag MS ( M ) 


u 






1.065(22) 


-0.034(16) 


0.794(21)(2) 


d 






-0.344(14) 


-0.034(16) 


-0.289(16)(1) 


s 


V 


= 32 3 64 





-0.031(12) 


-0.023(10)(1) 


T 3 


L p 


2.33 fm 


1.409(24) 





1.082(18)(2) 


r 8 






0.721(26) 


-0.006(18) 


0.550(24)(1) 


£ 






0.721(26) 


-0.098(42) 


0.482(38)(2) 


u 






1.071(15) 


-0.049(17) 


0.787(18)(2) 


d 






-0.369( 9) 


-0.049(17) 


-0.319(15)(1) 


s 


V 


= 40 3 64 





-0.027(12) 


-0.020(10)(1) 


T 3 


L - 


a 2.91 fm 


1.439(17) 





1.105(13)(2) 


T 8 






0.702(18) 


-0.044(19) 


0.507(20)(1) 


£ 






0.702(18) 


-0.124(44) 


0.448(37) (2) 



Nevertheless, having simulated only at one lattice spac- 
ing and sea quark mass, we cannot extrapolate our results 
to the ph ysic al point. Consequently, we underestimate 
the value [32| qa = 1.2670(35) from neutron /3-decays by 
13 % and find AT 3 = 1.105(13) instead. Our prediction 
AT$ = 0.507(20) differs by the same 13 % from the phe- 
nomenological estimate [32j as = 0.585(25). We take this 
as an indication of the size of the remaining systematics 
and add an additional 20 % error to all our results. 

Conclusions. — We determined the first moments of 
proton flavour singlet and non-singlet polarized parton 
distributions from nf = 2 lattice QCD, at a pion mass of 
285 MeV, at a single lattice spacing a « 0.073 fm. We 
found AS = Au + Ad + As = 0.45(4)(9) and a small 
negative As = -0.020(10) (4), in the MS scheme, at a 
scale \x = \/TA GeV. We underestimated both $a and 
a 8 by similar factors « 0.87 and this may suggest that 
some of the systematics cancel when considering ratios of 
matrix elements. Nevertheless, we emphasize that there 
is a considerable uncertainty in the as value (Io| and our 
AS is already relatively large, due to the small difference 
AT 8 - AS = -3As = 0.059(29)(12). 

Interestingly, our results are in remarkable agreement 
with the cloudy bag model prediction of 11] . The small 
(unrenormalized) As lat value obtained recently in [3l[ 
is also consistent with our study. Our AS value is 
larger than previously expected, however, it is compati- 
ble with the latest COMPASS number a ( % /3GeV) = 
0.35 (3) (5). The experimental number may increase fur- 
ther once smaller cc-values become accessible. We sug- 
gest relaxing the weak hyperon decay SU(3) f constraint 
on as in determinations of polarized parton distribution 
functions j and including our As prediction instead. 
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